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Abstract. We consider full perturbations to a covariantly defined Schwarzschild 
spacetime. By constructing complex quantities, we derive two decoupled, covariant and 
gauge-invariant, wave-like equations for spin-weighted scalars. These arise naturally 
from the Bianchi identities and comprise a covariant representation of the Bardeen- 
Press equations for scalars with spin- weight ±2. Furthermore, the covariant and gauge- 
invariant 1+1+2 formalism is employed, and consequently, the physical interpretation 
of the energy- momentum perturbations is transparent. They are written explicitly 
in terms of the energy-momentum specified on spacclikc three-slices. Ultimately, a 
Cauchy problem is constructed whereby, an initial three-slice may be perturbed by an 
energy-momentum source, which induces resultant gravitational fields. 
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1. Introduction 

In this paper, we consider perturbations to a background Schwarzschild spacetime which 
is described covariantly in terms of three non-vanishing locally rotationally symmetric 
(LRS) class II scalars [7j. Using the LRS class II symmetries, the background may be 
foliated with a family of spacelike three-slices [10] , and every three-slice is further foliated 
by a family of two-slices. The energy-momentum is specified as a first-order perturbation 
to an initial spacelike three-slice, and a Cauchy problem is constructed. This energy- 
momentum is what drives, or induces, any subsequent perturbed gravitational fields. A 
interesting example can be constructed by specifying a particular astrophysical object 
as the source to any resultant gravitational radiation. 

We use the 1+1+2 system of equations governing the first-order geometric and 
kinematic quantities which arise from both the Bianchi and Ricci identities. In this 
paper, we generalize the system given in [7] to include energy-momentum sources. 

We then demonstrate how to derive covaraint wave-like equations arising from the 
once-contracted Bianchi identities. Upon inspecting the fully non-linear equations, 
we note that the successful decoupling method described in [5] is also suitable for 
the system that arises here. Thus we construct complex tensors for the irreducible 
parts of the gravito-electromagentic (GEM) fields. Subsequently, by using a covariant 
decomposition of the two-sheet described in [5], we derive two decoupled wave-like 
equations for complex spin-weighted [H] scalars. These are a covariant representation 
of the Bardeen-Press (BP) equations pp. Furthermore, the energy-momentum source 
terms have been included to allow pragmatic calculations of perturbed gravitational 
fields arising from initial energy-momentum perturbations. 

2. Decomposing the Spacetime 

The necessary mathematical tools for decomposing tensors are separated into three 
distinct parts. The first represents the 1+3 formalism which is well established 
throughout the literature (for example, see [HI [HI E] ) - The second focuses on a further 
decomposition using the recently developed 1+1+2 covariant two-sheet formalism [7j. 
Finally, a further covariant decomposition of the two-sheet, expanding the work of [5], 
is described. As this paper uses the methods developed in [5], we retain that notation 
throughout this manuscript!]]. 

2.1. 1+3 Decomposition 

The tools required to irreducibly decompose the equations governing general relativity 
into a covariant 1+3 form are presented here. An observer is moving in the direction of 

\ The transformation between the notation used in this manuscript to that of [H El [10] is 
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V <llu> :=± V {flu) - - 1_^1_ Q/3 V Ql3 , (3) 



a timelike vector which is normalized such that 

g a/3 n a n p := -1. (1) 

The projection tensor is used to project any quantity onto the instantaneous rest space 
of the observer [9] and is therefore defined, 

-L^ := 9nu + n^n u . (2) 

Let V^ v represent any tensor, then the projected, symmetric and trace-free (PSTF) part 
of V^ v is defined using angular brackets, 

1 
3 

where _L represents an operator which projects all free indices, and the standard 
symmetric (round) brackets have been employed. A spacelike tensor, otherwise referred 
to as a three-tensor, will vanish when any of its indices are contracted with n M . 

The fundamental geometric quantities are found from the decomposition of the 
covariant derivative of the timelike vector, 

V„n„ = -A^ — ^ ± fU , K — n^h u + e^/up. (4) 

Here, K := — _L a/3 V a n^ is the expansion and := — V <il n v> is the shear. The 
spacelike acceleration is denoted n M := n a V a n M , and in general a dot derivative of 
any tensor is V^. M := n a V a V^, M . The spacelike vorticity is := \ a np, and 
the completely anti-symmetric three-Levi-Civita psuedo-tensor is defined such that 
tfiua '■= e a^ua na , where e^vor is the usual Levi-Civita pseudo tensor. 

The trace-free part of the Riemann tensor, otherwise know as the Weyl conformal 
curvature tensor (C Mi/cr7 -), can be expressed in terms of two spacelike, trace-free gravito- 
electromagnetic (GEM) tensors |10j . 

E^ u := C ai ,p u n a n p and := ~e afl/3y C^ ' wv n a n" '. (5) 

Here E^ v is the gravito-electric field and B^ u the gravito-magnetic field. 

Finally, the energy-momentum tensor (T^v) rnay also be irreducibly decomposed 
into 1+3 form, 

p := T a/3 n a n^ P := 1 l_ a/3 T aP , Jfl := - ± T^n a and 7T^ := T <fiu> . (6) 

Here p is the mass-energy density, P is the isotropic pressure, is the spacelike mass- 
energy flux and 71^ represents the spacelike, trace-free, anisotropic pressure. 

Before proceeding to the next decomposition, it is important to define two primary 
derivatives used throughout. The Lie derivative in the direction of a vector field is 
denoted Cx- Let any tensor be denoted V^,, M a "' T and therefore 

- \),.,- n : v n .\" - ... - r,,.,/" "v„.y 

+ V a ...S~ T V li X a + ... + V^. a °- T V„X a . (7) 
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In the case of a scalar (if) this reduces to a directional derivative 

£ x <p := X a V aV >. (8) 

It is noted that the Lie derivative of a covariant, spacelike tensor is also spacelike. 

The covariant three-derivative (-D M ) has the restriction that it operate on 
spacelike quantities only. Let W tl ... v <r "' T represent any three-tensor (W fl „.„ <T "' r n M = 
0, . . . , W^J* ' r ri T = 0), then the covariant three-derivative is defined 

n~w ■ ' -.=± v-u;,.,/'- ; . (9) 



2.2. 1+1+2 Decomposition 

The further decomposition of the 1+3 quantities into 1+1+2 form is now presented. This 
is a covariant, and gauge-invariant, decomposition recently developed by [7J, however, 
we again adhere to the notations and conventions employed in [5|§|. Another vector is 
defined (iV^) which is both spacelike (N a n a = 0) and normalized according to, 

+ a/3 N a N p = 1. (10) 

Therefore, another projection tensor that now projects all quantities (orthogonal to both 
and n/ 4 ) onto a two-sheet j7] is, 

S^:=±^-N^N U . (11) 

Furthermore, following the work of [7J, the PSTF part of some tensor (V^,), with respect 
to N a , is defined 

Vfav} '■= V&p) ~ -^S^ivS al3 V a j3, (12) 

where a "bar" over an index implies that the index is projected onto the two-sheet, i.e. 

V^:=S«S/V aP . (13) 

A tensor which vanishes when any index is contracted with both and is referred 
to a two-tensor. 

The important geometric quantities arise from the irreducible 1+1+2 decomposition 
of the covariant three-derivative of iV M , 

D„N V = -A„ v - -S^ K + N fJi N u + f e aP . (14) 

The expansion of the two-sheet is /C := —S a ^D a Np. The shear of the two-sheet is 
Afj, v '■= ~D{ fl N u y, which is symmetric and trace-free with respect to the two-sheet, i.e. 
S a/3 A a /3 = 0. We also denote := N a D a N^, and in general the "hat" derivative 
of any three-tensor (W /1 ...„ <r ™' r ) is defined W ^ v a - T := N a D a W^ u a - T . The scalar 
£ := ^e a/3 D a Np represents a twisting of the two-sheet, and the anti-symmetric two- 
Levi-Civita tensor is defined such that e^ u := e aflu N a . 

§ Where different, the transformation from the notation used in this pa- 
per to that of [7] is {/C, iV M , Sfj, v , K, E, S^, S^, A^ui B, 8^,8^, N^} -» 
{— 4>, n^, Npv, — 0, — E, — E M , — E^, — C/ii/j — W, — TCfj,, —Ti^y, a^}. 
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As expressed in [TJ, the 1+3 geometric and Weyl quantities are decomposed into 
the irreducible set {K, A, Q, E, £, B, A^, f^, E M , 8^, B^, E^, B^} according to 

A^ = cfy + ^4n M , (15) 
n M = ^ + ^AT M , (16) 
u^ = Q^ + QN^, (17) 

V = E^ - ~ S^E + 2 E {fl N v) + E AyV„, (18) 

^ = ^ - o 5 ^ + 2 + 5 ( 19 ^ 

= B»v-\ S^B + 2 B [fl N u) + B AyV (20) 

Here A := —n a N a = N a h a , a M := A^, := hp, fl^ := and Q := N a u a . 
Furthermore, E := A a/3 N a N^, E M := y4^ Q A^ a and E^ := A^j, and now E M and 
E Mi , are both two-tensors. Analogous definitions follow for both and B^. 
The 1+3 energy-momentum quantities, in 1 + 1+2 form, are 

J^J. + JN^ (21) 

7i> = - ^ s^n + 2 + n Ayv,, (22) 

where :— jp, J :— N a j a and the decomposition for TT^ is analogous to ({TBI . Finally, 
the covariant two-derivative associated with the two-sheet (dp) must operate on two- 
tensors only. Let U^, M a '" T represent any two-tensor, then 

d^J^- T := D^Up.J- f . (23) 
2.3. Decomposing the Two-Sheet 

This final section gives a fuller description of the covariant decomposition of the two- 
sheet as presented in [5j. The decomposition of the two-sheet which is suitable for 
decoupling the Bianchi identities follows a different trend as the 1+3 and 1+1+2 
decompositions. A natural decomposition of the two-sheet may be expressed by defining 
a complex-conjugate pair of vectors, m M and m M , which satisfy the following relationships 

m 



m a m a = 1, m a m a = 0, m a m a = 0. (24) 

This complex-conjugate pair are orthogonal to both n M and iV M . Consequently, they can 
be raised and lowered using S^ v . The projection tensor associated with the two-sheet 
may now be covariantly decomposed according to 

S lw = 2m( Jt ffi u ). (25) 

Consider the arbitrary two-tensors and U^, i.e. a contraction of any of their 
indices with n M or will vanish. By using ( |25i) . and recalling that U^ v can always 
be split into a symmetric plus an anti-symmetric part, they are irreducibly decomposed 
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according to 

Up = (U a m a )m^ + (U a rh a )m^ (26) 

+2m rt ([/ a/J m'V). (27) 

The generalization to tensors of more complicated type is clear. 

We now discuss the concept of spin- weighted objects [H]. Let the complex null 
vector, m M , undergo a transformation on the two-sheet according to 

m^CC^m^, (28) 

where C is an arbitrary complex scalar field and C its complex conjugate. Then any 
quantity, C^... u a '" T , which has a corresponding transformation of 

C,,.,/'--- -CX-'C,,...,/'--'. (29) 

is said to have a spin-weight s defined]]], 

s:=\{p-q). (30) 

Let us now reconsider the decomposition of according to ( 1261) . Since this is a two- 
tensor, there are only two independent scalars arising, namely (U a m a ) and (U a rh a ). 
Under the transformation fl28l). we have 



{U a m a ) ^CC- 1 {U a m a ), (31) 

thus (p, q) = (1,-1) and it therefore has a spin- weight of s = 1. An analogous 
consideration shows that (U a rh a ) has spin- weight s = — 1. Similarly, in the 
decomposition of U^, there are four independent scalars i.e. U a pm a m^ , U a pfh a fh° , 
U a pfn a mP and 6 r a/ gm Q m /3 , which have respective spin-weights of s = 2, s = —2, s = 
and s = 0. Thus one must be aware of how the transformation properties of spin-weight 
scalars differ from the usual zero spin-weighted scalars. 

We now consider the important quantities which arise from the decomposition of the 
covariant two-derivative of the complex-conjugate vectors. All the subsequent equations 
will naturally occur in conjugate pairs. However, we only display one of the pair and 
note that the other is found by taking the complex conjugate. Thus we have 

dpm v = -m ll m v rh' x Xa - ^ S^rrf'xa ~ i ^ & ( 32 ) 

where 

Xf, := m a d^rh a = ~Xn (33) 
is purely imaginary and has zero spin- weight. We also define 

a := im a fnpe al3 from which it follows a 2 = 1. (34) 
Thus a has zero spin- weight and since a = a, it is purely real. 
|| The boost weight (b) can also be defined, b — \(p + q)- 
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Therefore, (I3"2"j) depends only on the complex-conjugate pair (m a x a , m°Xa) which 
have spin-weights 1 and —1 respectively. We also have a constraint and a relationship 
for the divergence, which are respectively 

(d a + x a ) m a = and e a,3 d a mp = iad a m a . (35) 

Finally, using (1271) the two- Levi- Ci vita tensor is decomposed as 

e IM/ = i2am\ JJ ,rh u ]. (36) 

3. Background Schwarzschild Spacetime 

In the background spacetime, all the vorticity terms vanish, i.e. Q = and = 0. 
Therefore, under the conditions of vanishing vorticity only, the timelike vector n M is now 
normal to a spacelike three-slices and according to Frobenius theorem, this implies 
must satisfy 

n^VuTia] = 0. (37) 
The most general solution to (J3"T|) is 

n M = -QtV A »/, (38) 

where a is the lapse function and / is the Cauchy time function such that aV M / > 0. 
We now introduce two more vectors: t M which is tangent to the temporal coordinates 
world line and a spacelike shift vector /3 M :=_L The tangent is decomposed according 
to 

r = + (aC t f)n tl . (39) 

Therefore, the Lie derivative (in the directon of n M ) operating on a covariant and 
spacelike tensor becomes, 

£n= —^-t (A -£/?)■ (40) 

Similarly, the twisiting of the two-sheet in the background also vanishes i.e. £ = 0. 
Therefore, under the additional condition of zero twisting, the two-sheet meshes to form 
a two-surface, where is the corresponding normal. Since we require this normal to 
be hypersurface orthognal, we again use Frobenius theorem to conclude 

N, = 1 D )1 h, (41) 

where 7 and h are analogous to the lapse and Cauchy time functions such that 
^D^h > 0. A spacelike vector, r M , is defined to be tangent to the radial coordinates 
worldline and a shift vector for the two- surface is defined := r^. Then r M is 
decomposed according to 

r' l = U i +( r fC r h)N li . (42) 

Finally, the Lie derivative (in the direction of iV M ), when operating on a covariant two- 
tensor becomes 

C N = ^—{C r -C b ). (43) 
7 L r h 
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Hitherto, (I3"7)) - (|4"3"|) have only assumed that the background spacetime have vanishing 
vorticity and twisting. They constitute a covariant description of how the Lie derivatives 
and normal vectors are aligned with respect to coordinate tangents in vorticity/twist 
free spacetimes. 

We now impose further restrictions to arrive at a covariant description of the 
Schwarzschild spacetime. It was shown in [7J [11] that the background Schwarzschild 
spacetime may be covariantly described by only three non-vanishing LRS class II scalars 
(-4, /C, £), for which their governing equations are 

(C N + A-K)A = 0, {C N -A-^IC)IC = and S = AJC. (44) 
4. First-Order Perturbations 

We introduce a perturbation operator which is denoted 5. Subsequently, a perturbed 
tensor V^.„„ may be approximated to first-order according to, 

%...u « V»...v + SV^..,,, (45) 

where V^, M is the background tensor, and 8V^_ M is the first-order component of order 
e, where |e| << 1. 

It is important to stress that a background three-tensor will have a corresponding 
first-order tensor which is not spacelike with respect to the background. Instead, it will 
attain a first-order timelike component (2j HJ [13] . Similarly, a background two-tensor 
will have a corresponding first-order tensor which is not a two-tensor with respect to 
the background. 

Consider some perturbed two-tensor which therefore satisfies the covariant 
relationships U a N a = and U a n a = 0. These relationships, which also hold true in 
the background spacetime (U a N a = and U a n a = 0), are then expressed to first-order 
according to 

U a 5N a + N a 5U a = and U a 5n a + n a 5U a = 0. (46) 

Thus, it is clear that the first-order term (SU^) is not a two-tensor with respect to the 
background (i.e. does not vanish when contracted with n M or iV M ), and instead it now 
satisfies the first-order relationships ( )46l) . 

Now consider the perturbed two-tensor, U fJi ^ M . This may be expanded to first-order, 
and its corresponding first-order component is decomposed into 1+1+2 form in the usual 
fashion, 

Ufj,... v w U^.v + SU^ 

« U,...v + SU^p + (N a 5U a ... u ) N. + ... + (N a 5U^. a ) N v 

- {n a 5U a ... v )n^ - ... - {n a 5U^.. a )n u . (47) 

In the situation where the background two-tensor {U^, M ) vanishes, the first-order terms 
reduce to two-tensors according to (j4"6"j) . We then reuse the background symbol by 
replacing 5L^... P with U^ M , and similarly for scalars. 
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With the above definitions and conventions, the first-order geometric and Weyl 
quantities are, 

This set of first-order variables is gauge-invariant under infinitesimal coordinate 
transformations. This is a consequence of the Stewart- Walker-Sachs lemma [T5], [16], 
which states that a quantity which vanishes in the background has a corresponding 
first-order quantity that is gauge- invariant. There are three quantities which do not 
vanish in the background, and will experience first-order increments according to 

£^/C + 5/C, A&A + 5A and 8^8 + 58. (49) 

Thus <5/C, 5 A and 88 are not gauge- invariant. Finally, the first-order energy- momentum 
perturbations are all gauge-invariant, as they vanish on the background, hence 

{ P , p, j, j„ n, n M , iv) = o(e). (50) 

5. 1+1+2 Identities 

We now consider the Bianchi identities. It is the twice contracted Bianchi identities 
which provide the linear equations governing the energy-momentum perturbations. It is 
the once-contracted Bianchi identities for which the covariant, decoupled BP equations 
are naturally constructed with energy-momentum sources. Finally, the Ricci identities 
for both and provide the remaining covariant equations governing the full 
dynamics of the perturbed spacetime. 

5.1. Twice Contracted Bianchi Identities 

The twice contracted Bianchi indentities give rise to the equations governing the 
conservation of energy-momentum, V a T^ a = 0. They may be decomposed into 1+1+2 
form to give evolution and propagation equations^, which after linearization, become 

C nP + (C N + 2A-JC)J + d a J a = 0, (51) 

C n J + (C N + A - - K.)U + d a U a + (C N + A)P + Ap = 0, (52) 

CnJp + (C N ~ + d a Ii m - X - d^U + d„ P = 0. (53) 

These equations (along with additional equations of state) govern the first-order 1+1+2 
energy-momentum quantities. Furthermore, these equations clearly decouple from the 
remaining Bianchi and Ricci identities, as there are no first-order geometric quantities 
present. Thus, since the background quantities are prescribed, the system fl5Tl) - fl53l 
forms an initial-value problem and can be completely solved prior to solving for the 
remaining perturbations. These energy-momentum perturbations therefore become 
known sources in the following identities. 

% dEE]) derives from 7^V Q T MQ = 0, ([52]) from iV^V Q T MQ = and ([53]) from V"l)- iQ = 0. 
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5.2. Once Contracted Bianchi Identities 

The once-contracted Bianchi identities for Einstein's equations have a rich 
correspondence with the equations governing electromagnetism. We show how two 
decoupled, covariant and complex equations arise naturally from these identities. 
Therefore, we give them extra attention and provide a discussion of the fully non-linear 
equations before writing down the linearized system. The identities can be expressed in 
terms of both the Weyl conformal curvature and energy-momentum tensor (T^), and 
we define 

B vtTT := V a C avaT - {V [a T T]v + ^g v[a W T] T) = 0. (54) 

The 1+3 equations governing the GEM tensors can be found by an irreducible 
decomposition of (1311) . which ultimately gives rise to two constraint and two evolution 
equations^], 



D a Efia + CuapB^AyP + 3 B^ a LO a 

= 1 ZV - ± D a K m + 1 Kj^ - X - A«j a + ^ e^jccjp, (55) 



D a B IMX — e iUX pE ori Arf — 3 



= (p + P) uo^ + ^e^DJp + ~ e^n^A^ - ~ n„ a cu a , (56) 

(£„ - ^ K)E <flv> + e af3 ^(D a + 2 h a )B„f + 5 A <fl a E u>a - e { ^E u)a u p 
= ~ (p + P) A^ - - {D <(1 + 2 n <fi )j v> 

- l - {C n + ~ K)7T^ U> - 1 7T a<tl A v> a + \e { ^ v)a u p (57) 

and 

(Ai - I K)B <fiV> - e aP{fM (D a + 2 ti*)E v f + 5 A <iX a B v>a - e { ^B v)a u p 

= X ~e a ^A v ff -^e a ^D a 7V v f - ^j<^ v> . (58) 

In the absence of energy-momentum sources, it is easy to inspect that these fully 
non- linear 3+1 Bianchi identities ()55|) - (|58p are invariant under the simultaneous 
transformation E^ v — > B^ and B^ u — > —E^ v . Therefore, following the work of [5], a 
natural way to decouple the system is to construct a complex-conjugate pair of the form 
E^u ± iB^ u . It is then straightforward to reintroduce the energy-momentum quantities 
into the calculations. 

We now proceed to write the linearized Bianchi identities in a 1+1+2 form. 
Therefore, we define three first-order 1+1+2 complex tensors, 

^■^Sepp + iSBps,, Qp^SEp + iSBp and 6$:=6£ + i6B, (59) 

+ The constraint (|55|) is derived from _L n a n^ B a p^ = or _L £[ M Q ^ B^ a p = and (|56|) from 
e^ 1 n a B a p 1 = or _L n a Br Att/ i , = 0. The evolution equations, ([57]) and ([58]) . derive from 

_L n a B^ v ) a = and _L £( M a ^ B v ) a p = respectively. 
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where and $ M are gauge-invariant, whereas is not because SS because in the 
background $ = £ ^ 0. Furthermore, when constructing these complex quantities, it is 
found that some of the other geometric quantities also naturally combine, and thus we 
make three additional definitions 

A M „ := <5£ A p + i e a (t ,6A p)Q} A„ := 5N p + i (e M ° <5£ a + Stop), (60) 

T^:=^ + ie%M a - (61) 

With these definitions, (|55l)-( l58|) . under a first-order perturbation, and in 1+1+2 
form areEl, 

6[(£ N - -)C)$}+d a ® a -3ien = g, (62) 
(C N - /C)$ A + - ^ 5(d^) + ^SA, = g„ (63) 

5(£ n $) - i e^fy + ^ (E - 2 -K) - i 2 i = T, (64) 
£ n $ A + ^ £ T M - i e ( /cF$ 7)Q - i ~ [<f(d a $) - (2 .A + /C)$ a ] = ^(65) 

+ i e^/^ + \k, + 2 A)$ p)a - i e { /d w $„ } - -£ A^ = (66) 

where the vertical bars in ( l66l) imply that the index is excluded in the symmetric 
brackets. The first-order energy-momentum, or forcing, terms have been defined 

G := \ C N {p - ^n) + ^/Cn-i d a U a + i ~ e at3 d a J p , (67) 
Q„ ■= I d^p + - rf M n - 2 ( £ at - /C)n A - - d a Ii m - i - e«(C N Ja - (68) 

.F : = -I £„n - \ (C N + 2 A + I K) J + \d a J a - i 1 e a ? djip, (69) 
2 o 2 6 2 

^ := -~£J^-^(£ iV + 2A+2/C)J A -~d Ml 7 

-i ^ e M a Q rf a p + 2 d a U - (C N - X - £)n a J , (70) 



1 1 .1 



(71) 



Each individual term in the 1+1+2 Bianchi identities fl62l) - fl66l) has been carefully 
grouped such that they remain gauge-invariant under infinitesimal coordinate 
transformations. For example, in (164"|) the first term, <5(£„$), is gauge-invariant because 
its corresponding background value vanishes, i.e. £ n $ = 0. 

The energy-momentum sources (|671)- (17T1) appear as large complicated terms. 
However, the dynamics of all the energy-momentum terms can be fully prescribed using 
the conservation of energy- momentum equations ( l5TT) -(|53l). 

* These equations are derived as follows: ([62]) from (7V 7 n Q n' 3 + ie' 37 rt Q; )i? a ^ 7 =0, ([63)) from n a rfi 'B + 
ie^n a B af3l = 0. PI) from N N^n a B Pia + i N"< e af) B^p = 0. ([65]) from (N v n a B {fiv)a + 
i N v ^B v)afi = and © from n Q S (pp)Q + i e^B^ = 0. 
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5.3. 1+1+2 Ricci Identities 

The remaining 1+1+2 equations arise from the irreducible decomposition of the Ricci 
identities for both n M and iV M , which are defined respectively 

Qv» a ■= 2 V^V^n* - R^arn T = 0, (72) 
R^, a ■= 2V [M V H iV CT - i^iNT = 0. (73) 

It is also useful to note that the Riemann tensor may always be expressed as a 
sum of its trace (R^ u and R) and trace-free parts (C^t). Then, by substituting the 
Einstein equations to eliminate both the Ricci tensor and Ricci scalar in terms of the 
energy-momentum tensor, it becomes 

2 

We now present the first-order system and note that the following equations ( 1751) - 
( |94|) generalize those presented in [7J to include energy-momentum terms, for which 
we separate these into four groups. 

5.3.1. Constraints for the two-slices These are derived as follows: (1751) derives from 
N^rfR^ = 0, (|75|) from e^n a R^ ua = and ([77]) from S ua R-^ a = 0. 

- rf M (S - ^K) + 2 e/t^fl - 2 cPE^ + /C S M + K e„ a n a - 2 e«B a = X - (75) 
d a tl a -e af3 d a £ p - (2A + )C)n+B = 0, (76) 

-i 5(<yc) - e/<U + d Q A- + ^ = n M . (77) 



5.3.2. Constraints for the three-slices These are derived as follows: f[78|) derives 

from N»S UT Q^ a = 0, (J7JD from e^Q^ = 0, (JH0D from N^R^ = 0, ([HID fr o m 

S^N^R^v = 0, (1521) from n^n a R^ a = 0. (1531 from N^R^y = and @3D from 
iV^ (W ) = 0. 

- | AT) - | /C £ + - e^tkfy, = J7, (78) 
(£tv - -4 - /C)fi + (TO, = 0, (79) 
(^-/C)£-^d a iV> = 0, (80) 

- „4 - i/C)/C] - - X/C) + d Q iV Q = ~ p + I n, (81) 
/^(^ + e A a fi a ) - /C S M + ePE^ - l - d M (£ + ^tf) - efd a n + 2 < Ae M a fi Q = 0, (82) 

£n-A^ + d^N u y — S^ u = - 11^, (83) 

(C N + ix:)E^ - d {M (E„> - e v} a n a ) - e ( /^ )Q = 0. (84) 



Covariant Schwarzschild perturbations I 



13 



5. 3. 3. Transportation equations These are derived as follows: (1851) from _L M<T Q/j, ua n u = 



0, dSOD from ep^RafrnP = and (JHTD from n a N^R aPfi = 0. 

C n K + 6[(C N + A- K)A] + d a A a = 1 (p + 3 P) , (85) 

e/£ n O Q - ^ £jvA a + X - 5{d^A) - X -A N„ = 0, (86) 

-£ n A> p + (£* + A)a A - (A+ ~/C)(E„ - e/fi Q ) - e, a 8 a = ~ J,. (87) 

5.3-4- Evolution equations These are derived as follows: n^N^Q ^ Da ^ = 0, (I8"9l 

from n^^Q^ = 0, © from r^i^} = 0, (JSTJ) from n^Q^ = 0, © from 
^5 W V = 0, (JHSD from n^e ua R^ va = and (f5¥|) from n^N y n a R^ va = 0. 

£„(S A + e/^) + 5(d„A) + (.A + ^ /C)A M = — n M , (88) 

£ n Q-^e Q ^ Q ^-Ae = 0, (89) 

£nA- - + g /C ) S ^ + + = °> ( 90 ) 

t~"tv^-'j±v ^/if A A^ -|- (ij^Ajy} — n^i/, (91) 

5(£ n /C) + (A + ^ /C) (S - H at) + ^ = (92) 

- (A + \ JC)tl - ~ e a0 d a a p + l& = 0, (93) 

£n(S - ^K) - d a A a -5{£-AJC) = -~(p + 3P)-~Il. (94) 



For the coming decoupling section, we require evolution and propagation equations 
for the complex tensor A^. For the propagation equation we simply take the appropriate 
combination of (183]) and (1841) . Similarly for the evolution equation we use both (190]) and 
( 19T]) . They are respectively 

— £jvA^p — i e^ a — - K + i e^ a d v yA a = i - e^ a Il u ) a , (95) 

1 1 

-£ n A^p + -iie^A^ - i-/Ce (M Q S jy)a + ie { ^ a d u} T a = -11^.(96) 

6. Decoupling the Bianchi Identities 

With the once contracted Bianchi identities expressed in complex form, it is now possible 
to derive a decoupled equation for Before, proceeding with the derivation, it 

is necessary to first present the commutation relationships for the three derivative 
operators when acting on any covariant first-order tensor vV--^> 

£ndp,Pv...a ~ d fl C n ^ u ,„ a = 0, (97) 
£>Ndp,ipv...cT - d^C N ^p D ...„ = 0, (98) 
(C N + A)£ n (Pii...p - C n C N Lp n ,„ D = 0. (99) 
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The calculation is complicated, lengthy and begins by taking the Lie derivative 
with respect to of fl66l) . Then commute the two Lie derivatives, which are acting 
on and substitute ( 1661) again to eliminate £ n $^p. Also commute C n d^ acting on 
$ M and substitute (|65|) to eliminate £ n $^. Then similarly commute C^d^ and then 
substitute (1631) to eliminate £tv$ /j . Throughout this entire procedure, collate any terms 
which have a background gravito-electric scalar (£). Finally, this collated term must 
be reduced using both fl95l) and fl96l) . This results in a covariant, and gauge-invariant, 
wave-like equation 

[£„£„ - (C N + A-K, )C N - d a d a - V) $ AP + i (4 A + 2 /C)e ( /£ n $ p)Q = M^. (100) 
The potential and source have beed defined respectively, 

V := -2K g ~^IC 2 -4A 2 + 2S, (101) 
3 

Mfj, v := CnT^y + i e^ a (C N — A — - /C)J^) Q - i e^ a d\ a \ T v } + rf {A1 ^}. (102) 

Here, K g is the Gaussian curvature of the background two-surface, i.e. K g := | ^R, 
where &>R is the Ricci scalar for the background two-surface. 

We now proceed to further decompose the wave-like equation for <&^ u f llOOl) into 
two decoupled wave-like equations for spin- weighted scalar s. Both $^ and M.^ are 
symmetric and trace- free and they are decomposed using ( l27l) . 

= X rh fl rh u + Y m^m u , (103) 
= Qm^ffiv + Pm^my. (104) 

Here, the scalars X and P have spin- weight s = 2, whereas F and P have spin- weight 
s = —2. Thus, from the decomposition of (11001) . we get two decoupled, covariant and 
gauge-invariant, wave-like equations for the spin-weighted scalars X and Y, 

(£„ + 4 .4 + 2 /C)£ n X - (£iv + A - 3 /C)£jvX - (d a + 4 x a )d Q X 

- [-2K 9 + /C 2 -4^ 2 + 4 X ]X = g, (105) 

and 

(£„ - 4 4 - 2 /C)£ n F -(C N + A-3 K)C N Y - (d a - 4 x Q )d Q y 

- [-2K 9 + /C 2 -4^ 2 + 4x]F = P, (106) 

where x '■= {d a wP){d a rhp) = % These wave-like equations correspond to the BP 
equations p] for gravitational perturbations^. These covariant and gauge-invariant 
wave-like equations may be solved for the spin-weighted scalars (X and Y) and then 
reconstructed back into <&^ u which is zero weighted. Finally, the real and imaginary 
parts give respectively 6£^ u and 5B^ V . 

jt These equations have been checked with the vacuum equations in the literature [IJ [6l [12] , by using 
a specific choice of the definitions for the null and complex conjugate- vectors as defined in [6]. 
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The purpose of this section is to express the wave-like equations ( 11051) and ( 1106}) in 
a pragmatic coordinate form suitable for both analytic and numerical computation. 
The covariant equations governing the Schwarzschild spacetime (jHJ) are solved for the 
background LRS class II scalars [7] to give 



„ M , 
A = ~2 I 3 



2M 



2 / 2M 
K = — \1 

TV V 



and 



£ 



2M 



(107) 



This corresponds to the standard static Schwarzschild line element (ds) using 
coordinates x M := (t, r, 9, 0), 



ds z 



1 



2M 



dt 2 + 



-dr 2 + r 2 (d6 2 + sin 2 6d(j) 2 ) , (108) 



where M is the Schwazschild mass. Without loss of generality the Cauchy time function 
is set equal to t and the lapse function is therefore, 



a 



1 - 



2M 



(109) 



Furthermore, the complex-conjugate vectors defining the two-metric (<S M „) is chosen such 
that 

' (O,O,l,icosec0), (110) 



V2i 



and the Gaussian curvature of the two-spheres and other necessary quantities reduce to 
1 . , cot 2 6 



K, 



a=l, 



Xn = [0, 0, 0, i cos 9} and 



X 



(111) 



The shift- vectors {j3^ and 6 M ) vanish and the Lie derivatives ( 1401 and ( j43l) when operating 
on scalars reduce to 



r 

*~"n 



1 d 



and 



Cm = a 



d 



a dt dr 
Therefore, the wave-like equations fll05|) - (11061) become 

1 d 2 fl2M-Ar\ d 



(112) 



a 2 dt 2 



2Mr dt 



X 



+ 



a 



d 2 (10M-Qr)d 2 fr 2 — QMr + QM 



and 



dr 2 
1 d 







d 2 



r 2 - 2Mr 
4 cos 6* d 



X 



~^ 2 di 2 + 



sin 2 9 d(p 2 
12M-4r\ d 



+ i 



sin 2 9 



+ 4 cot 2 9 



X 



-Q 



(113) 



- 2 -2Mr <9t 



F 



+ 



cu 



<9 2 (10M-Qr)d 2('r 2 -QMr + QM 
shi^^r 111 00 J + sin 2 90 2 1 sin 2 9 dd> 



r 2 - 2Mr 
. 4 cos 9 d 



Y 



+ 4 cot 2 9 



Y = -P. 



(114) 



Covariant Schwarzschild perturbations I 



16 



These wave-like equations (I113p and fjl 14|) give the exciting prospect of calculating 
the first-order gravitational fields resulting from a first-order energy-momentum 
perturbation. They constitute an initial- value problem, whereby the energy-momentum 
is specified as a perturbation to an initial spacelike three-slice as the source to resultant 
gravitational fields. In particular, they are suitable for modelling gravitational radiation. 

An initial hypersurface may be constructed at t — to an d the energy-momentum 
perturbations are specified 

{ j0j p 1 j I j,,n ) n, 1 n (1 „}| t=l0 . (115) 

The conservation of energy- momentum equations (I5T]) - (I53]) . plus equations of state, may 
then be used to evolve the energy-momentum perturbations for all time. In the case of 
numerical computations, they can be integrated up until the events of interest. 

Thus, with the energy-momentum sources (Q and P) now known, we further require 
initial conditions for X and Y, i.e. 

dY 



x\ dx 



and Y 

t=t 



t=t , m 



(116) 

t=t 



Since we know m M and m M , this corresponds to initial Cauchy data given by 

(117) 



<9$ 

3Vlt=i and 



dt 



t=t 

To obtain this initial data, one must solve the non-trivial first-order initial-value 
problem. That is to say, one must solve the first-order system arising from the Bianchi 
and Ricci identities (!75|) - (!M|) on the initial spacelike three-slice. It is also important to 
separate the particular solutions arising from the energy-momentum perturbations from 
the homogeneous solutions which exist irrespectively. 

Finally, in the absence of energy-momentum sources (P — Q — 0), the homogenous 
solutions to (11 13)) and (1 1 14j) encompass the usual vacuum perturbations, whereby the 
vacuum operators are separable (as there is a common factor of 1/r 2 out the front of all 
angular terms) and spin-weighted spherical harmonic solutions arise 



7. Summary and Discussion 

We have shown how constructing complex quantities for the GEM tensors was 
successful for finding decoupled, gauge- invariant and covariant, wave-like equations 
(11131) and (II 14j) . These equations form an initial- value problem and are suitable 
for calculating first-order gravitational perturbations which are induced by first-order 
energy-momentum sources. 

The Weyl conformal curvature tensor comprises ten trace-free components of the 
Riemann tensor. The two spin-weighted complex scalars X and Y correspond to four of 
the ten components. That is, X and Y correspond to through (11031) which further 
corresponds to £^ LU and through ([59]) . Now since the two-tensors and are 
trace-free with respect to the two-surface, they each have two components for a total of 
four. 
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The study of the remaining six components is the subject of future research. 
Immediate difficulties arise when one attempts to derive wave-like equations for the 
remaining complex GEM tensors. Consider $ At , by taking the Lie derivative with respect 
to of f|65|) . one is confronted with a term involving £ n T^, which does not arise 
naturally. Finally, 5$ is not gauge-invariant and thus the physical interpretation of this 
scalar is nebulous. 
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